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Soft value calculation for multilevel signals 

This invention relates to digital communications systems and, more 
particularly, the generation of soft reliability values for multilevel signals. 

5 

Within the field of digital communications, multilevel modulation is used to 
map a number of bit sequences to a signal alphabet comprising a number of 
signal symbols, i.e. a number of points in signal space. For example, a bit 
sequence may be mapped onto a point in a complex signal space. A signal 

10 alphabet of size M allows log2(M) bits to be mapped to each symbol. 
However, when symbols are received at a receiver, they may be affected by 
noise, thereby affecting the decoding of the signal when retrieving the 
transmitted bit sequence. If multilevel modulation is used in conjunction with 
channel coding, many channel decoders, such as iterative decoders based 

15 on the BCJR algorithm, require so-called soft bit values as an input. A soft bit 
value corresponds to a reliability value of a single bit being 0 or 1 . 

Examples of multilevel modulation include multi-amplitude level modulation in 
Pulse Amplitude Modulation (PAM), multi phase level modulation in Phase 
20 Shift Keying (PSK), multi signal point modulation in Quadrature Amplitude 
Modulation (QAM). 

For example, an emerging technology for wideband digital radio 
communications of Internet, multimedia, video and other capacity-demanding 

25 applications in connection with the third generation of mobile telephone 
systems is the evolving Wideband Code Division Multiple Access (WCDMA) 
specified as part of the 3GPP standardisation organisation. Within this 
technology, High Speed Downlink Packet Access (HSDPA) is provided 
including a high speed downlink shared channel (HS-DSCH) which uses 

30 16QAM. In 16QAM for example, M=16, i.e. each symbol in the signal 
alphabet represents 4 bits. Future releases may comprise even larger 
constellation sizes such as 64 QAM. 

It is known how to convert signal symbols to soft bit values by calculating all 
35 distances in signal space between the received symbol and all signal points 
of the signal alphabet. In particular, in order to obtain optimal performance, a 
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likelihood ratio is calculated depending on corresponding sums of 
probabilities where the probabilities are functions of the calculated distances. 
It is further known that in the calculation of a likelihood ratio the sums of 
probabilities may be approximated by the dominant contributions to the sums 
5 of probabilities in the likelihood ratio (A.J. Viterbi, "An intuitive justification and 
a simplified implementation of the MAP decoder for conventional codes", 
IEEE Journal on selected areas in communications, 16(2), February 1998). 

Even though this approximation significantly reduces the computational 
10 complexity while only causing a negligible loss in performance, a calculation 
of all distances to all the signal points is still required in order to determine 
which two are actually needed for the calculation of the likelihood ratio. For 
example, in a 16QAM modulation, 16 distances have to be calculated for 
each received symbol. In particular, if a high rate of symbols needs to be 
15 decoded, e.g. several hundred symbols per millisecond, the above 
performance issue is particularly severe. 

German patent no. DE 199 12 825 describes a method of receiving a data 
symbol in which the closest constellation point of a corresponding 8-PSK 

20 multi-symbol constellation is identified for a received data symbol. 
Subsequently, for each bit position, the symbol having a different value at 
that position and which is closest to the identified constellation point is looked 
up from a look-up table. A soft value is calculated as the difference of the 
distances from the received symbol to the two identified constellation points. 

25 Hence, the number of distance calculations is reduced. 

It is an object of the invention to provide a method of generating a reliability 
value that reduces the computational complexity. 

30 The above and other problems are solved when a method of generating a 
reliability value for a received multilevel signal in relation to a number of 
predetermined signal symbols each associated with a corresponding bit 
sequence including a first bit position; the reliability value being indicative of 
likelihood information of receiving said multilevel signal comprises the steps 

35 of 
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- identifying a first one of the number of signal symbols as being closest to 
the received multilevel signal; and 

- estimating the reliability value based on a stored pre-computed distance 
function of at least the first signal symbol and a second one of the number 

5 of signal symbols, where the second signal symbol is the signal symbol 
closest to the first signal symbol of the signal symbols corresponding to a 
different binary value at the first bit position of the respective associated 
bit sequence than the first signal symbol. 

10 Given a signal constellation, .for each signal symbol and for each bit position 
it is known which one of the other signal symbols having an opposite value at 
that bit position is closest to that signal symbol. Based on this information, a 
distance function of the signal symbol and the closest signal symbol with 
opposite bit position may be pre-calculated and stored in a look-up table. 

15 Here, the term pre-computed distance function comprises any intermediate 
result of the computation of the reliability value from the first and second 
signal symbol. 

It has been realised by the inventors that the calculation of the reliability 
20 values may be based on a pre-computed distance function where the 
function is pre-computed given the first signal symbol and a second signal 
symbol, i.e. a pre-computed function of the first and the second signal 
symbol. The second signal symbol is the signal symbol that is closest to the 
first signal symbol of all signal symbols of the number of signal symbols that 
25 have a different binary value at the first bit position of the respective 
associated bit sequence than the first signal symbol. Hence, the complexity 
of the calculations to be performed upon receipt of a signal symbol is 
considerably reduced. Hence, the above method is particularly well-suited for 
low-complexity implementations in mobile receivers, as it reduces the 
30 required computational resources. 

It is a further advantage that the closest signal point needed for calculating 
the likelihood ratio is determined first, such that only the corresponding 
distances to the identified signal points need to be determined. In this way, 
35 the computational complexity is reduced significantly, since not more than 
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two distances need to be calculated in order to determine the likelihood ratio 
for each bit. 

According to a preferred embodiment of the invention, the stored pre- 
5 computed distance function comprises the distance between the first signal 
symbol and the second signal symbol, and the step of estimating the 
reliability value further comprises the step of determining a first distance 
between the received signal and the first signal symbol. Hence, the actual 
distance between the first and the second symbol is pre-calculated and 
10 stored. Once the first signal symbol is identified, this stored distance may be 
looked up and used as an approximation for the distance between the 
received signal and the second signal symbol. Consequently, a further 
reduction in computational complexity is achieved, since only one distance 
has to be calculated. 

15 

When the step of estimating the reliability value comprises the step of 
determining a polynomial function of the first distance and the second 
distance between the first signal symbol and the second signal symbol, 
multiplied by a predetermined constant, the computational complexity is 
20 further reduced, as no logarithm needs to be calculated. In one embodiment, 
the polynomial function is a difference of the squared distances. 

According to yet another preferred embodiment of the invention, the stored 
pre-computed distance function is indicative of one of a number of functional 

25 relations between the received multilevel signal and the reliability value, and 
the step of estimating the reliability value further comprises the step of 
selecting a functional relation of said number of functional relations 
dependant on the first signal symbol and the first bit position. Hence, the 
calculation of the likelihood value only comprises the step of calculating the 

30 corresponding stored function. Preferably, the functional relationship is a 
linear function of a signal component, thereby reducing the calculation to a 
multiplication operation and an adding operation. 

According to yet another preferred embodiment, the stored pre-computed 
35 distance function comprises, for each signal symbol and bit position, an 
approximation of the corresponding reliability value. Hence, according to this 
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embodiment, an approximation of the reliability value may be directly looked 
up once the closest signal symbol is identified, thereby providing a 
computationally very efficient method which eliminates the need of online 
distance calculations. 

5 

Preferably, the stored information is stored in a look-up table comprising a 
plurality of pre-computed distance functions indexed by the number of signal 
symbols and the bit positions of the number of bit sequences, thereby 
providing fast access to the information. 

10 

In one embodiment of the invention, the method further comprises the step of 
providing the reliability value as an input to a decoder, e.g. an iterative 
decoder using the BCJR algorithm or any other decoder using soft values as 
an input It is an advantage of the invention that it provides an accurate and 
1 5 resource-efficient approximation of soft values as an input to such decoders. 

The first signal symbol may be identified by comparing the signal 
components with predetermined thresholds or decision boundaries, for 
instance by means of a slicer, i.e. a circuit which compares a signal with 
20 predetermined thresholds. Hence, a fast and computationally inexpensive 
method is provided for identifying the closest signal symbol without the 
necessity of calculating all distances between the received value and all 
signal symbols. 

25 It is a further advantage of the invention that cost-effective, standard 
components may be employed when implementing a method according to 
the invention. 

When the likelihood information comprises a log-likelihood ratio, a high 
30 performance quality is achieved, as the use of a log-likelihood corresponds a 
theoretically optimal way of calculating soft reliability values. However, other 
methods of calculating likelihood information may be employed, such as a 
log-likelihood of the signal power. 



35 In a preferred embodiment of the invention the step of identifying the first 
signal symbol as being closest to the received multilevel signal comprises the 
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step of identifying the first signal symbol as being closest to the received 
multilevel signal with respect to a Euclidean distance measure in a signal 
space, as the Euclidean distance is directly related to the probabilities of a 
likelihood calculation. Alternatively, other suitable known metrics may be 
5 used instead Euclidean distances. 

The signal space may be a real or complex signal space. For example, in 
QAM modulation two amplitude-modulated signals are transmitted on a 
single carrier, but shifted in phase by 90 degrees. Hence, the resulting signal 
10 points may be represented in the complex plane representing the so-called 
in-phase (I) and quadrature (Q) components of the QAM signal. In general, in 
M-QAM, the corresponding signal constellation comprises M signal symbols, 
where M=2 n , n=2, 3, 4, 5, 6, 7, etc. 

15 When the number of signal symbols is associated with the number of bit 
sequences such that the bit sequences associated with all nearest 
neighbours of each signal symbol only differ from the bit sequence of that 
signal symbol at one bit position, the error rate of the transmission system is 
reduced. This form of mapping is referred to as Gray mapping. 

20 

The invention further relates to an arrangement for generating a reliability 
value for a received multilevel signal in relation to a number of predetermined 
signal symbols each associated with a corresponding bit sequence including 
a first bit position; the reliability value being indicative of likelihood information 
25 of receiving said multilevel signal; characterised in that the arrangement 
comprises 

- first processing means adapted to identify a first one of the number of 
signal symbols as being closest to the received multilevel signal; 

- storage means adapted to store information related to the first signal 
30 symbol and a second one of the number of signal symbols being closest 

to the first signal symbol and corresponding to a different binary value at 
the first bit position of the respective associated bit sequence than the first 
signal symbol; and 

- second processing means adapted to estimate the reliability value on the 
35 basis of the stored information. 
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The arrangement may be implemented by any processing unit, e.g. a 
programmable microprocessor, an application-specific integrated circuit, or 
another integrated circuit, a smart card, or the like. The term processing 
means comprises general- or special-purpose programmable 
5 microprocessors, Digital Signal Processors (DSP), Application Specific 
Integrated Circuits (ASIC), Programmable Logic Arrays (PLA), Field 
Programmable Gate Arrays (FPGA), etc., or a combination thereof. The 
processing means may be a CPU of a computer, a microprocessor, a smart 
card, a SIM card, or the like. The first and second processing means may be 
10 separate processing means, e.g. separate circuits, or they may be combined 
in one processing means, e.g. performed by suitable instructions executed by 
a programmable microprocessor. 

The term storage means includes magnetic tape, optical disc, digital video 
disk (DVD), compact disc (CD or CD-ROM), mini-disc, hard disk, floppy disk, 
15 ferro-electric memory, electrically erasable programmable read only memory 
(EEPROM), flash memory, EPROM, read only memory (ROM), static random 
access memory (SRAM), dynamic random access memory (DRAM), 
synchronous dynamic random access memory (SDRAM), ferromagnetic 
memory, optical storage, charge coupled devices, smart cards, etc. 

20 

Furthermore, the above discussed features and steps of the method 
according to the invention may be incorporated in the above arrangement 
according to the invention. 

25 The invention further relates to a device for receiving multilevel signals 
comprising an arrangement as described above and in the following. 

The device may be any electronic equipment or part of such electronic 
equipment, where the term electronic equipment includes computers, such 
30 as stationary and portable PCs, stationary and portable radio 
communications equipment. The term portable radio communications 
equipment includes mobile radio terminals such as mobile telephones, 
pagers, communicators, e.g. electronic organisers, smart phones, PDAs, or 
the like. 

35 
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The invention will be explained more fully below in connection with preferred 
embodiments and with reference to the drawings, in which: 

figs. 1a-b schematically show a receiver according to an embodiment of the 
5 invention; 

fig. 2 shows an example of a signal constellation with 16 signal symbols; . 
fig. 3 shows a flow diagram of a method of determining reliability values; 

10 

fig. 4 shows an example of a look-up table for use in the method of fig. 3; 

fig. 5 shows a flow diagram of a method according to an embodiment of the 
invention; 

15 

fig. 6 shows an example of a look-up table according to the embodiment of 
fig. 5; 

fig. 7 shows another example of a signal constellation with 16 signal 
20 symbols; 

fig. 8 shows a flow diagram of a method according to another embodiment of 
the invention; and 

25 fig. 9 shows an example of a look-up table according to the embodiment of 
fig. 8. 

Fig. 1a schematically shows a receiver according to an embodiment of the 
invention receiving a radio signal from a transmitter via a communications 

30 channel. Transmitter 101 is adapted to send a signal s via a noisy channel 
102 to receiver 103. The signal s represents one of a set of M signal points 
Si...S M in a signal space where each signal point is related to a respective bit 
sequence of log 2 (M) bits. In the presence of noise in the transmission 
channel 102, the receiver 103 receives a signal f that deviates from the 

35 transmitted signal s. In one embodiment the signal is a Code Division 
Multiple Access (CDMA) signal using a spread spectrum technique. The 
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receiver 103 comprises a receiver circuit 107 for transforming the received 
spread spectrum signal into the signal symbol r. The receiver further 
comprises a channel decoder 106 for decoding the received signal symbol r, 
e.g. a BCJR or Viterbi decoder. The decoder 106 requires soft bit values as 
5 an input. Hence, the receiver 103 further comprises a circuit 104 which is 
adapted to calculate soft values for the log2(M) bits of the received signal 
symbol r and to provide the calculated soft values to the decoder 106. 
According to the invention, the receiver 103 further comprises a memory 105, 
such as on-chip memory, EPROM, flash memory, or the like, in which a look- 
10 up table is stored for use in an efficient calculation of the soft values by the 
circuit 104, preferably as described in connection with figs. 3-9. 

Fig. 1b schematically shows a more detailed block diagram of the receiver 
circuit 103 of fig. The circuit 107 comprises a RAKE receiver 110 suitable for 

15 receiving CDMA signals, i.e. a receiver which uses several baseband 
correlators to individually process several signal multipath components. The 
correlator outputs are combined to achieve improved communications 
reliability and performance (see e.g. "Digital Communications" 4th Edition, by 
John G. Proakis, McGraw-Hill, 2000). The sampled received radio signal r* is 

20 fed to the RAKE receiver 110 which generates the signal symbol r to be 
decoded. The circuit 107 further comprises a channel estimator 111 and a 
noise estimator 112, e.g. implementing any suitable channel estimation and 
noise estimation technique known in the art. The channel estimator receives 
the received radio signal r\ identifies up to N different radio paths or channel 

25 taps and estimates corresponding delays Ak, k=1 N, and complex channel 

estimate h f = (hri,...h,N) of these paths. The channel estimator 111 further 
provides a set of complex combiner weights w = [wi w 2 .. . w N ] T to be used by 
the rake receiver. Here, D T denotes a transposed vector. For example, the 
weights may be determined according to an optimisation criterion, such as 

30 maximising the received signal energy. 

The calculated delays A k and the combiner weights are provided to the rake 
receiver 110. The RAKE receiver 110 comprises delay circuits 115 which 
delay the incoming signal according to the N channel taps. Further, the 
35 receiver 1 1 0 comprises circuitry 1 1 6 for multiplying the N delayed versions of 
the received signal with a spreading code c for dispreading the spread 
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spectrum signals and circuitry 117 for summing the signals to form a radio 
symbol. Furthermore, the rake receiver 110 comprises multiplier circuitry 118 
for multiplying each of the N radio symbols with the combiner weights (w k )*, 
k=1,...,N, where ()* denotes complex conjugation. Finally, the RAKE receiver 
5 110 comprises an adding circuit 119 which combines the weighted symbols 
to form the received symbol estimate r which is fed to the soft value 
calculation circuit 104. 

When using a multilevel constellation of signal points Si...Sm the amplitude 
10 information should be maintained in order to ensure successful demodulation 
in the receiver. Consequently, the reference points Si...S M should be scaled 
properly. In the following it is assumed that the channel estimator 111 
estimates the channel gain on the basis of a reference channel h r which has 
a channel gain that may be different from the actual gain of the traffic 
15 channel, e.g. a HS-DSCH. The gain difference between the reference 
channel and the traffic channel may be denoted with g. Hence, the received 
symbol r after the RAKE receiver 110 may be expressed as 

r=g w H h r s + n, 

20 

where w H is the Hermitian conjugate of w, w H h r denotes an inner product, s 
is the transmitted symbol and n is a noise term, e.g. representing additive 
white Gaussian noise (AWGN). The gain parameter g is signalled to the 
receiver, w is selected by the combiner in the receiver, and h r are the 
25 channel estimates. Hence, at the receiver, the reference signal symbols 
Si...S M may be scaled appropriately, according to 

Sj^gw^Sj, j = 1,...,M. (0) 

30 In fig. 1b, the receiver circuit 107 comprises circuit 113 adapted to calculate 
the above scaling factor g w H h r and a multiplier circuit 114 for multiplying the 
reference symbols with the scaling factor, resulting in properly scaled signal 
symbols S 1f ...,S M , which are fed into the soft value calculation circuit 104. 

35 Finally, the noise estimator 112 provides an estimate of the signal noise level 
a which is fed into the soft value calculation circuit 104. 
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According to the invention, the soft value calculation circuit locates the signal 
symbol which is closest to the received signal r and calculates corresponding 
soft values L m for bit m f m=1,...,Iog 2 (M), e.g. according to one of the 
5 embodiments discussed in connection with figures 3-9. 

It is noted that the receiver circuit described in connection with figs. 1a-b 
merely serves as an example, and the scope of the invention is not limited to 
the type of receiver, nor to the above scaling of signal symbols. 

10 

Fig. 2 shows a signal constellation with 16 signal symbols. The signal 
constellation comprises M=16 signal points Si through Si 6 in a two- 
dimensional signal space, e.g. the l/Q components in a 16QAM signal 
constellation. Preferably, the signal points are distributed regularly, such that 

15 the distance to the nearest neighbours of each signal point is the same. The 
reference points may take values that suit the implementation in question. 
However, alternatively, other signal constellations may be chosen. In fig. 2, 
16 different bit sequences 0000 through 1111, each consisting of log 2 (M)=4 
bits, are mapped onto the signal points S r Si 6 . Preferably, the mapping of the 

20 bit sequences to the signal points is chosen such that the bit sequence of 
each signal point only differs from those of the nearest neighbours by one bit, 
thereby optimising the decoding performance. For example, in fig. 2, signal 
point Ss has three nearest neighbours, S 4 , S 7 , and S12. The bit sequence of 
S 4 , i.e. 1111, differs from the sequence 1110 of S 8 only at bit position 4, etc. 

25 Alternatively, other mappings may be chosen. 

For every bit position m mapped on a signal point, the signal points in the 
constellation may be divided into two sets where the signal points in each set 
have bit value 0 and 1, respectively, at that position. In the following, the set 
30 of signal points with a 0 at the m-th position is denoted Ao, m , and the 
corresponding set with a 1 at the m-th position is denoted Ai, m . For example, 
in the example of fig. 2 and for m=1, A 1t1 ={S3, S 4 , S 7 , S 8 , Sn, Si 2 , S15, Si 6 } 
and Ao,i={Si, S 2 , S 5l S 6 , S 9 , Si 0> S 13 , S i4 }. The sets are of equal size with 
M/2 elements each. 

35 
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When one of the signals Si,...,Si 6 is transmitted over a noisy channel, the 
received signal will differ from the transmitted signal according to a 
corresponding distribution. The actual shape and width of the distribution of 
received signals depends on the characteristics of the noise. In fig. 2, the 
5 cross 201 represents an example of a received signal r. 

Prior to providing the received 16QAM radio symbols to a decoder, e.g. a 
turbo decoder, they are converted into soft values. Hence, a soft value is 
calculated for each bit of every 16QAM symbol. A soft value of the m-th bit in 
10 the sequence mapped to r may be defined as 

P(r|Sm=l) (1) 



P(s m =l)K r 


l s m=l) 


P(s m =0)P(r 


|s m =0) 



15 where s m is the m-th bit in the bit sequence represented by the transmitted 
signal, and P(s m =i| r), i=0,1, are the a posteriori probabilities of the bit s m 
where r is the received signal. It is noted that the second equality assumes 
that s m =1 and s m =0 are equally probable in the chosen alphabet. Otherwise, 
the overall ratio of probabilities should be taken into consideration in the 

20 following. However, this would only give rise to a constant factor. Hence, L m 
corresponds to a log-likelihood ratio of probabilities. The probabilities 
P(r|s m =i) in eqn. (1) may be written as 

P(r |s m =i)=^ £P(r,s) i = 0,1. (2) 
25 S€A im 

Hence, the calculation of the above probability involves a summation over 
M/2 terms each including a joint probability P(r,s). This is a computationally 
expensive task, especially if M is large, e.g. M=64. 

30 In many applications, the above soft values L m may be approximated by 

max P(r,s) (3) 

I Inn S£A ^ | nn P ( r l§1,m) 

SsA f 



35 
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where §i, mi i=0,1, are the signal points that result in the largest contribution to 
the sums in eqn. (2). Hence, in the calculation of the probabilities, the sums 
over M/2 terms are approximated by the their respective dominant terms, 
according to 

5 

log ]TP(r,s)«log max P(r,s) (4) 

The above approximation is often referred to as the "max log MAP" 
approximation which yields a good approximation in cases where the above 
10 sums are dominated by one term, as for example in the case of Gaussian 
noise when the signal to noise ratio (SNR) is large. The above probabilities 
depend on the distances between the received signal r and the respective 
signal points. For example, in the case of additive zero-mean Gaussian noise 
with variance o 2 , the log-likelihood ratio of eqn. (3) may be expanded as 



15 



a"" 2 



exp(-|r-§ 1|m | 2 /a 2 ) . 

L m-log } (l r - § 0.m| -Ml.m| ) (5) 

a- 2 exp^|r-§o,mr/^ 2 J V ; 

20 It is noted that it is assumed that the si, m are scaled corresponding to the 
received signal according to equation (0) above. In the following, we define 
di,m=|r- Si, m |, for 1=0,1, to be the distances between the received signal r and 
the closest signal points in the sets Aj, m , respectively. For example, in fig. 2, 
for m=1, d 0l i corresponds to the distance 82 between r and the closest signal 

25 point in Ao.i, i.e. S6, while du corresponds to the distance 81 between r and 
the closest signal point in Ai f i, i.e. Ss. According to the invention, the 
likelihood ratio in equation (5) is obtained by first identifying the closest 
signal point S 8 , and then determining the distances 81 and 82, as will be 
described in greater detail below. Hence, a computationally expensive 

30 calculation of all the distances between r and all the signal points S1...S16 in 
order to identify the shortest distances 81 and 8 2 is avoided. Alternatively to 
identifying S 6 by means of a look-up table and then calculating 82, the 
distance 83 between S 6 and S 8 is looked up and used as an approximation 
instead of 8 2 , thereby saving additional computational resources. This will be 

35 described in greater detail in connection with figs. 5-6. A further embodiment 
of the invention will be described in connection with figs. 8-9. 
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It is noted that, preferably, in the above estimation of the reliability values, a 
proper scaling of the signal points in the QAM constellation is taken into 
consideration. If this scaling is taken into consideration, the above log- 
5 likelihood ratio may be written as 

L m =a- 2 {|r- §0>m | 2 -|r-§ 1>m | 2 } <«> 



10 As described in connection with fig. 1b, when using a multilevel constellation 
as in the example of fig. 2, the amplitude information should be maintained in 
order to ensure successful demodulation in the receiver. Consequently, the 

reference points Sj, j=1 16 should be scaled properly. If this scaling is 

taken into consideration, the above log-likelihood ratio may be written as 



15 



25 



U^K^If-So^l 2 -!?-^! 2 ] (7) 



20 r = 



i.e. with the properly scaled signals 
r 



gw"h r 



^=-%r m=1 i°g 2 (M),/=o,i, (8) 

and where 

K = (gw H hJ/a 2 (9) 
is a constant which depends on the signal to noise ratio. 



30 Fig. 3 shows a flow diagram of a method of determining reliability values. 
According to this method, the soft values L m are calculated using the 
approximation in equation (7). Initially, in step 301 , a signal r is received and, 
in step 302, the signal point s from the set of signal points Si...S M which, 
according to a Euclidean metric, is closest to r is identified. For example, an 

35 efficient way of identifying s is by means of a slicer. In step 303, the distance 
5i between r and s is calculated. Subsequently, for bit positions 
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m=1 log 2 (M), the following steps are performed: In step 304, the signal 

point § which is closest to s is looked up in a look-up table 308. This signal 
point corresponds to the signal point which is closest to r and has the 
opposite bit value at position m than s. In step 305, the distance 62 between r 
5 and s is calculated. Based on the distances 5i and 82, the soft value L m is 
now approximated according to eqn. (7) above: If the bit value s m of s at 
position m is 0, the soft value is approximated by L m = K(52) 2 -(Si) 2 (step 

306) . Otherwise the soft value is approximated by L m = K (Si) 2 -^) 2 (step 

307) . Here, K is a constant which depends on the noise distribution as 
10 described above. Referring to the example illustrated in fig. 2, the closest 

signal point to the received signal r (marked by the cross 201) is s=S 8 . When 
calculating a soft value Li for the first bit position m=1 using the method of 
fig. 3, the first bit in S 8 is identified to be s<i = 1. From a pre-computed look-up 
table, e.g. as illustrated in fig. 4, the closest signal point with a "0" in the first 
15 bit position is s=S 6 . Hence, the distances 81 and S 2 may be calculated as 
Si=|r-S 8 | and S 2 =|r-S 6 |, respectively, where 1 1 denotes the Euclidean distance. 
Thus, the soft value Li is approximated by Li = K (5 2 ) 2 -(8i) 2 . 

Consequently, the method described above requires at the most 1 + log 2 (M) 
20 distance calculations, since the closest distance 81 has to be calculated once 
for a received signal (step 303) and, for each bit position m, the distance 8 2 is 
calculated in step 305. This is to be compared with M distance calculations 
when all distances to all signal points are calculated. Hence, the 
computational complexity of this method only grows logarithmically with the 
25 size of the symbol alphabet rather than proportional to the alphabet size. This 
is a considerable reduction of the computational complexity, in particular for 
large alphabet sizes. 

Fig. 4 shows an example of a look-up table for use with the method of fig. 3. 

30 The look-up table 308 identifies the pre-computed closest signal points and is 
indexed by the bit numbers m and the signal points Si...Sm. Each row 
corresponds to one of the signal points Si...S M . For example, row 402 
corresponds to signal point S 2 , such that each element in row 402 identifies a 
signal point which is the closest to S 2 among all signal points having a bit 

35 value opposite to S 2 at bit position m. Each entry in table 308 consumes 
log 2 (M) bits for identifying one out of M signal points. Furthermore, the table 
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consists of M rows and log 2 (M) columns. Consequently, the table requires M 
[iog 2 (M)] 2 bits. For example for M=8 the memory consumption is 72 bits and 
for M=16 the memory consumption is 256 bits. 

5 Fig. 5 shows a flow diagram of a method according to an embodiment of the 
invention. Again, this embodiment utilises the approximation of equation (7) 
for the calculation of the soft values L m . As in the method of fig. 3, in the initial 
step 501, a signal r is received and, in step 502, the signal point s from the 
set of signal points Si...S M , which is closest to r is identified, e.g. by means 

10 of a slicer. In step 503, the distance 81 between r and s is calculated. 
Subsequently, for bit positions m=1,.. M log 2 (M), the following steps are 
performed: In step 504, the distance 8 3 between s and the signal point s 
which is closest to § and has the opposite bit value at position m is looked up 
in a look-up table 508. Subsequently, this distance 8 3 is used as an 

15 approximation for the distance 5 2 between r and s when approximating the 
soft value L m according to eqn. (7) above. Hence, if the bit value § m of s at 
position m is 0, the soft value is approximated by Lm= K (8 3 ) 2 -(8i) 2 (step 

506) . Otherwise the soft value Is approximated by L m = K (8i) 2 -(8 3 ) 2 (step 

507) . Again, K is a constant which depends on the noise distribution. 
20 Referring again to the example illustrated in fig. 2, the closest signal point to 

the received signal r is s=S 8 . When calculating a soft value Li for the first bit 
position m=1 using the method of fig. 5, the first bit in S 8 is identified to be §1 
= 1. From a pre-computed look-up table, e.g. as illustrated in fig. 6, the 
distance to the closest signal point with a "0 M in the first bit position is di. 8 =8 3 . 
25 Hence, the distance 81 is calculated as 8i=|r-S 8 | and 8 2 is approximated by 8 3 . 
Thus, the soft value Li is approximated by Li = K (53) 2 -(5i)2- 

Consequently, as the pre-computed distance S 3 is used as an approximation 
for 8 2 , the method according to this embodiment requires only one distance 
30 calculation, i.e. the calculation of 81 (step 503). Again, this is to be compared 
with M distance calculations when all distances to all signal points are 
calculated. 

Hence, it is 'an advantage of this embodiment that the computational 
35 complexity does not grow with the size of the symbol alphabet, thereby 
yielding a computationally efficient method of approximating soft values. 
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Hence, it is an advantage that storing the look-up table only requires little 
storage space. 

It is a further advantage of the method according to the invention, that it 
5 yields a good approximation of the soft values, thereby providing a good 
decoding performance. 

Fig. 6 shows an example of a look-up table according to the embodiment of 
fig. 5. The look-up table 508 comprises the pre-computed distances between 

10 the signal points in a constellation of size M. A distance d m ,k in table 508 
denotes the Euclidean distance between signal point S k and the closest 
signal point with opposite bit value at position m. Assuming that each 
distance is stored with a resolution requiring t\ bits, each entry in the table 
508 requires j\ bits. Furthermore, the table consists of M rows and log 2 (M) 

15 columns. Consequently, the table requires i\ M log 2 (M) bits. For example, 
for M=8 the memory consumption is 24i\ bits and for M=16 the memory 
consumption is 64tj bits. Hence, it is a further advantage of this embodiment 
that it requires little storage capacity. In an embodiment where the resolution 
of the pre-computed distances is higher than Iog 2 (M) bits, i.e. •n>log2(M), 

20 processing time is traded for memory space in comparison with the method 
of figs. 3-4. 

It is noted that additional storage space may be saved by only storing each 
distance once, i.e. in case the same distance appears in two or more entries 
25 of the table, a reference to that distance may be stored in one of the entries, 
instead. 

Alternatively, other layouts of a look-up table may be used. For example, in 
one embodiment, the look-up table may comprise all M(M~1)/2 mutual 
30 distances between the signal points Sr, thus requiring t|M(M-1)/2 bits of 
storage. However, for M>4 this embodiment requires larger storage capacity 
than the embodiment of fig. 6. 

Fig. 7 shows another example of a signal constellation with 16 signal 
35 symbols. As in fig. 2, the signal constellation comprises M=16 signal points 
Si through Si 6 in a two-dimensional signal space, e.g. the l/Q components in 
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a 16QAM signal constellation. The signal points are distributed regularly, 
such that the distance to the nearest neighbours of each signal point is the 
same. In the example of fig. 7, they are assumed to be selected such that 

5 S k = Xk + j yk, where x k , yk e [-3d, -d, d, 3d], k=1,...,M, 

where d is an arbitrary constant and where j 2 = -1. For example, d may be 
chosen to d=1. However, alternatively, other signal constellations may be 
chosen. 

10 

In fig. 7, 16 different bit sequences 0000 through 1111, each consisting of 

log 2 (16)=4 bits, are mapped onto the signal points Si S 16 . Preferably, the 

mapping of the bit sequences to the signal points is chosen to be a Gray 
mapping, i.e. such that the bit sequence of each signal point only differs from 
15 those of the nearest neighbours by one bit, thereby optimising the decoding 
performance. 

As above, for every bit position m mapped on a signal point, the signal points 
in the constellation may be divided into two sets Ao.m and A^, where the 
20 signal points in each set have bit value 0 and 1 , respectively, at that position. 

Fig. 8 shows a flow diagram of a method according to another embodiment of 
the invention. Again, this embodiment utilises the approximation of equation 
(7) for the calculation of the soft values L m . As in the embodiment of fig. 3, in 

25 the initial step 801 , a signal r is received. The received signal may be written 
as r = Re(r) + j lm(r) and, in the following the magnitude of the I- and Q 
components of r will be denoted by a = | Re(r) | and b = | lm(r) |, respectively. 
After the received symbol is combined in the combiner, in step 802, the 
signal point § from the set of signal points Si...Sm, which is closest to r is 

30 identified. In this embodiment it is assumed that the constellation of signal 
points corresponds to the constellation of fig. 7. In fig. 7, each signal point 
corresponds to a decision region where the decision regions are separated 
by a set of decision boundaries 701 through 706. Hence, the closest signal 
point s may be found by performing two comparisons of the inphase 

35 component and the quadrature component, respectively. For example, if 
Re(r) < 0 (decision boundary 705) and Re(r) <-2d (decision boundary 706) 
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and if lm(r)>0 (decision boundary 702) and lm(r)>2d (decision boundary 
701), the received signal lies in the decision region corresponding to Si, i.e. 
s = Si is the closest signal point. Subsequently, for each bit positions 
m=1,...,log 2 (16)=1,...,4, the soft value L m may be calculated using the 
5 approximation of eqn. (7), assuming proper scaling. Consequently, in this 
example the soft value Li for the first bit is 

LKSO = K(|-a+j b - (d+3j d)| 2 - |-a+J b - (-3d+3j d)| 2 ) = K (8ad-8d 2 ), 

10 Hence, in the above equation, instead of computing two distances squared, 
each involving a calculation of the type |x+jy| 2 , the soft value may be 
calculated by scaling the inphase amplitude a of the received symbol with 
8dK and, subsequently, by adding a constant -8Kd 2 . It is further noted that 
the constant d may be chosen as any suitable positive real number. 

15 

The remaining three soft values for a received signal in the decision region 
corresponding to Si are accordingly: 

L 2 (Si) = K(|-a+j b - (-3d+3j d)| 2 - |-a+j b - (-3d-j d)| 2 ) = K (-8bd+8d 2 ) 
20 U(Si) = K(|-a+j b - (-d+3j d)| 2 - |-a+j b - (-3d+3j d)| 2 ) = K (4ad-8d 2 ) 
USi) = K(|-a+j b - (-3d+3j d)| 2 - |-a+j b - (-3d+j d)| 2 ) = K (4bd-8d 2 ), 

as in the constellation of fig. 7 the closest symbols with opposite second, 
third, and fourth bit compared to Si are S 9 = -3d - jd, S 2 = -d + 3jd, and 
25 S 5 = -3d + jd, respectively. 

The table 808 of fig. 9 illustrates the calculated soft values for all decision 

regions corresponding to the symbols Si Si 6 , and for all bits, m=1,...,4. As 

can be seen from table 808, all soft values may be calculated by scaling one 

30 of the inphase component a or quadtrature component b of the received 
signal rand subsequently adding a constant. Hence, using the pre-calculated 
equations of table 808, the soft values may be calculated in a very efficient 
way. In one implementation, each entry of the look-up table 808 may 
comprise the scaling factor, the constant to be added and a bit indicating 

35 whether it is the inphase component a or the quadtrature component b of the 
received signal r which is to be scaled for a given soft value. Preferably, the 
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table is indexed by the decision region and the bit values. It is noted, 
however, that many of the entries of table 808 are identical. Consequently, it 
will be apparent to a skilled person that table 808 may be stored in a memory 
efficient manner, e.g. by storing a list of the distinct entries and, in table 808, 
5 referring to the corresponding list members. In general, it is noted that 
constellations which are Gray coded or show another regularity, the 
redundancy of the entries in table 808 may be utilised to reduce the memory 
consumption of table 808. 

10 Referring again to fig. 8, in steps 804-805 the soft values for the identified 
decision region and for all bits are calculated. In step 804, the relation to be 
calculated, i.e. the scaling factor and the constant to be added, are retrieved 
from a stored table 808 in memory, e.g. a look-up table as shown in fig. 9. 
The retrieved relation is calculated in step 805 resulting in the soft value for 

1 5 the corresponding bit number. 

It is noted that the processing load in the receiver may further be decreased 
by pre-calculating the relations of table 808 and by storing the pre-calculated 
soft values: Assuming that the inphase and the quadrature components of 

20 the received signal each are quantised to n bits, the soft values of table 808 
may be precalculated and tabulated for every different inphase and 
quadrature value, thereby further decreasing the required calculations, as the 
scaling and adding of step 805 are not necessary in this embodiment. 
However, such a table of pre-calculated soft values increases the memory 

25 consumption. Above, a=|Re(r)| and b=|lm(r)| were defined as the absolute 
values of the real and imaginary parts of r, respectively, i.e. without sign 
information. Hence, a and b, each are represented by n-1 bits. If each of the 
pre-calculated soft values is to be represented by m bits, the total memory 
consumption of a full table is m 2"" 1 4 16 bits (for each of the 16 decision 

30 regions and each of the 4 bits, 2°" 1 different soft values are stored, each with 
a precision of m bits). For example, for n=4, the total memory is 512m bits. 
Note, however, that in this embodiment, the pre-calculated table still needs to 
be multiplied with the factor K. 
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It is noted, that the above memory consumption may be further reduced by 
utilising the fact that many of the entries of table 808 are identical and by 
utilising the symmetry of the constellation of fig. 7. 

5 It is noted that the invention was described in connection with soft values 
defined as a log-likelihood ratio indicating a reliability value for the bit values 
of a received sequence. However, other definitions of soft values depending 
on the distance of the received signal to the signal points may be used as 
well. 

10 

It is further noted that the signal constellations of figs. 2 and 7 are merely 
used as examples. The calculation of soft values according to the invention is 
not limited to these signal constellations. 

15 Finally, it is noted that the embodiment described in connection with figs. 5-6 
is particularly well suited for large signal constellations, as it saves memory, 
whereas the embodiment of figs. 8-9 is particularly well suited for 
implementing a medium-size signal constellation, e.g., 16QAM, as it saves 
computational resources. 

20 



